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Abstract — Optimal precoding in the relay is investigated to maximize 
ergodic capacity of a multiple antenna relay channel. The source and 
the relay nodes are equipped with multiple antennas and the destination 
with a single antenna. It is assumed that the channel covariance matrices 
of the relay's receive and transmit channels are available to the relay, 
and optimal precoding at the relay is investigated. It is shown that the 
optimal transmission from the relay should be conducted in the direction 
of the eigenvectors of the transmit-channel covariance matrix. Then, we 
derive the necessary and sufficient conditions under which the relay 
transmission only from the strongest eigenvector achieves capacity; this 
method is called Maximum Eigenmode Relaying (MER). 



I. Introduction and Related Work 

Cooperative communication has been a rather active field of 
research in the recent past (e.g. ^ was ^ lst shown by van der 

Meulen in |5] that cooperation can enhance the transmission rate of 
a communication system. Later on, substantial work was carried out, 
investigating the effect of cooperation in various types of communica- 
tion systems. The first works concentrated on communication nodes 
with single antennas and several relaying protocols were proposed 
of which two important ones are (non-regenerative) AF (Amplify- 
and-Forward) relaying and (regenerative) DF (Decode-and-Forward) 
relaying. In this paper, we will focus on non-regenerative relaying 
techniques. 

Multiple antenna systems are well known to boost the Shannon 
capacity of a communication system (see e.g. j6)). A natural setup is 
the combination of cooperative communications with multiple anten- 
nas, which is at the heart of current research in the field. Assuming 
multiple antennas at the relay, one of the major tasks is to design a 
suitable relaying protocol. However, depending on the system or the 
desired performance criterion, different "optimal" relaying protocols 
can exist: for instance, the non-regenerative relaying protocols, e.g. in 
(7J-j9), are designed to minimize Mean Square Error (MSE) but 
other relaying protocols, e.g. in 1 1 1 — 1 13], are assumed to maximize 
Shannon capacity. Moreover, depending on the available channel 
information in the relay, the "optimal" relaying protocol can be 
different. 



A. Related Work 

A noncoherent cooperative system was investigated in |12| . The 
source and the relay nodes are equipped with multiple antennas and 
the destination node with single antenna. It was assumed that the 
source and the relay have access to the covariance matrix of the 
channels. Another assumption was that the antennas in the source 
are correlated but no correlation was assumed in the relay. Using 
that knowledge, the optimal transmit direction in the source and the 
relay was derived. Assume Q is transmit covariance matrix in the 
source with spectral decomposition Q — UqAqUq and S is the 
correlation matrix of the source, with spectral decomposition £ = 
U^A-sU^ . It was proved that the optimal Q must be in the form 
of Q° = UsAqoU^; i.e. Uqo = J7e and Aqo is in descending 
order which needs to be solved numerically. It was also proved that 
the optimal gain matrix in the relay is a weighted identity matrix 



which depends on ur, £ and Q. Furthermore, the necessary and 
sufficient condition for the optimality of beamforming in the source 
was derived as 

(I+7AFMI+7A?) 



7^2 < 



- 1 



(1) 



where V(l+y\f) is a function given in 

E {(1 + jz)e 1+! T(0, z)} where E{-} represents expectation opera- 
tion. Note that the evaluation of the optimality of beamforming in 
{I} requires computationally expensive monte carlo simulations or 
numerical integrations due to AfjyXf). 

More recent results were obtained in [13 ] where an optimal relay 
precoding was investigated for a system with correlated antennas at 
the relay. It was assumed that the relay has access to full CSI of Hi 
but only to the covariance matrix of H2, i.e. to Rr with spectral 
decomposition .Rr = UrArUr . In j 1 3| , |l4) , it was proved that 
the optimal relay precoding matrix, F°, is in the form of F° = 
UrA° f 'V , ^ 1 , where V h x is the unitary matrix with columns as the 
eigenvectors of HiH^ . Along with the numerical methods to derive 
Ap in 1 13 1 , the optimality of beamforming was considered only for 
the asymptotic case of high transmit SNR. The relay beamforming 
was found to be optimal if following inequality holds: 



1) < 1 



(2) 



a" o 
Details can be found in |13[ Sec. III-C]. 

In this paper, we consider a MIMO cooperative system with 
multiple antennas at the relay in which the relay has access to the 
covariance matrices of preceding and following channels. Moreover, 
the relay antennas are assumed to be correlated. With partial channel 
knowledge (the covariance matrices) available to the relay, a capacity 
maximizing relaying method will be introduced; the exact system 
model and the desired performance criterion will be discussed in 
forthcoming sections. 

The paper is organized as follows: In Section|n] the system model 
is explained and two problems are stated that are the main subject of 
this paper. Section [nij deals with designing optimal precoding in the 
relay in order to maximise capacity. In Section [IV] the necessary and 
sufficient conditions for the optimality of MER are investigated. In 
Section [V] some asymptotic results and their usefulness in practical 
systems are considered. In Section [VI] numerical results are presented 
and finally, some conclusion remarks are explained in Section [VII] 

II. System Model and Problem Statement 
A. Notation 

Matrices are represented by boldface upper cases (H). Column 
and row vectors are denoted by boldface lower cases (h), and 
hi indicates the i-th element of h. The superscript H stands for 
Hermitian transposition. We refer to the identity matrix by J. The 
expectation operation is indicated by E(-) and fx(x) is reserved for 
probability density functions (pdf); As represents a diagonal matrix 



2 



with elements organized in descending order and \f denotes the 
i-th diagonal element of As- For simplicity of notation, (\f) 2 is 
abbreviated by \f 2 . The trace of a matrix is denoted by Tr(-). 

B. System Model 

A dual hop, half duplex non-coherent MIMO communication 
system is considered in this paper. A source node with n$ antennas 
communicates with a single-antenna destination node only via a relay 
node that is equipped with n R antennas (that are used for both 
reception and transmission). It is assumed that a direct link between 
the source and the destination is not available. The half duplex 
constraint is accomplished by time sharing between the source and 
the relay; i.e. each transmission period is divided into two time slots: 
the source transmits during the first time slot and the relay during the 
second one. The relay remains silent during the source transmission 
and vice versa. It is assumed that the source does not have access to 
any statistical or deterministic channel state information (CSI). The 
signal received at the relay (y R ) due to the source transmission is 
given by 

y R = Hix + w R (3) 

where the tir x ns matrix Hi represents channel between the source 
and the relay (below, only the statistics of Hi are assumed to be 
known at the relay). With Ps the power constraint of the source, 
the column vector x is the signal transmitted from the source with 
Q = E(xx H ) — 0-fn s and the column vector wr represents the 
receiver noise in the relay with elements independently drawn from 
a complex Gaussian random variable with variance No. The relay 
multiplies y R with gain matrix F and forwards it to the destination. 
Then, the received signal at the destination is 



y D = h 2 Fy R + w D 

= h 2 FHix + h 2 Fw R + w D 



(4) 



where the row vector h 2 indicates the channel between the relay and 
the destination; wd represents the receiver noise at the destination. 
For simplicity, we assume that wo is statistically equivalent to wr 
and that both noise processes have unit- variance, i.e., No = 1; the 
latter choice is no extra restriction, as the ratios of transmit powers 
and noise powers determine performance, and we are still free to 
choose Ps and Pr arbitrarily. 

We assume spatial correlation only at the relay, which can be due 
to unobstructed relay node or space limits at the relay which force 
antennas to be closely located. Justifications to assume transceivers 
with spatial correlation can be found in (15) , |16|. The correlation 
matrix in the relay is represented by S with spectral decomposition 
£ = C7sAeC/e , where Us is a unitary matrix with its columns 
the eigenvectors corresponding to S, and As is diagonal matrix with 
the eigenvalues of £ in decreasing order. Note that we assume i.i.d. 
channels for Hi and h 2 . Therefore, the channel matrices Hi and 
h 2 can be written using the Kronecker model as 



Hi — S 2 i?i u 

i 

h 2 — h 2w T, 2 



(5) 
(6) 



where Hi w and h 2w are i.i.d., zero mean, unit variance complex 
Gaussian random variables, independent of each other; denotes 
a matrix formed by the square roots of the elements in the matrix S. 

C. Problem Statement 

With |4), the ergodic capacity of the system is defined as 
1 



C = 



max 
F-H\\FyR\\ 2 }<PR 



E{C(Hi,h 2 ,F)} 



(7) 



where the relay gain matrix F is to be designed to maximise 
when the expectation operation is carried out over Hi and h 2 . 
Assuming Q = !^-In s (this means equal transmit power from each 
antenna is chosen in the source, because no channel knowledge 
is available there), the channel capacity C(Hi, h 2 , F ) for given 
channel matrices is 



C(Hi,h 2 ,F) = log (1 



Ps h 2 FHiHi F H h 2 
ns N (l + h 2 FF H h?) 



(8) 



where No(l + h 2 FF H h 2 ) is the total equivalent noise power which 
remains constant per coherence time due to a block fading assumption 
we impose. 

Two major problems are addressed below: 

Problem I: Solving the optimization problem in {7}, in order to 
find the optimal precoder F° in the relay that maximizes mutual 
information between the transmit signal from the source and the 
received signal at the destination, given that the correlation matrix 
£ as well as the unit variance of elements of i.i.d Hi w and h 2m are 
available at the relay. 

Problem 2: Obtaining necessary and sufficient conditions based 
on £, so that maxium eigenmode relaying (MER) is optimal. 

III. Optimal Precoding in the Relay 
In this section, Problem 1 from Section|II-C|is addressed. However, 



before we continue to define F°, a closer look at the power constraint 
of the relay in |7| is provided. The power constraint of the relay is 
given by 

E{||Fy R || 2 } = E{Tr(Fy R y^F H )} < P R . (9) 

One can write Tr(Fy R y R F H ) = Tr(F H Fy R y R ). Note that F is 
a positive gain matrix; therefore, G = F F is a positive symmetric 
matrix. One can also choose the gain matrix F to be symmetric (i.e. 
F — G 3 ), which follows from substituting jij and jfijl in Q where 
F is left and right multiplied by S?; therefore, by choosing F = 
G 3 , the problem of finding F° will be replaced with finding G°. By 
substituting ([3j and <[5j in <[9j and applying spectral decomposition 
of S, l|9j will be further simplified to 

E{Tr(Gy K yg)} = P s Tr{GU s A s C/f ) + iV Tr(G) < P R . (10) 

By combining l | IQf and {7}, we have 
1 

- max 

2 ~ Ps r 



C : 



E{C(Hi,h 2 ,G)}. 



(11) 



G:P s Tr(G[/ E A E E/£ )+iV Tr(G)<PR 

Assuming |5]( and |6j, applying spectral decomposition according 
to E = J7eAe?7e and considering that the statistics of Hi and 
h 2 do not change by a multiplication with a unitary matrix, the 
"instantaneous" capacity ([8} for given channels Hi , h 2 is obtained 
as 



C(H 1 ,h 2 ,G) =log(l + 
where 7 = P&/Nq and 



jh 2w G 2 H i w Hi w G 2 h 2 
n s (l + h 2w Gh 2w ) 



') (12) 



G 2 = A 2 U"G^U S A 2 . (13) 
Then, the relay power constraint in ( |10) can be written as 

PsTr(A^ 1 G) + iV Tr(As 2 G) < Pr. (14) 



On the other hand, applying the spectral decomposition G = 
U qA-qU^, in |l2b and considering that the statistics of a random 
matrix do not change by multiplying with a unitary matrix, so we 
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obtain 



<7(JI ll h fll A <5 )=Iog(l + 



7 fe2 U ,A|g lm iJf m A|/if m , 



By careful inspecting \12\ and < | 1 5 [ > , the conclusion is that 
E(C(H 1 ,h 2 , G)) = E(C(iTi, fo 2 , A 6 )) . 



(15) 



(16) 



One interpretation of 1I61 is that C{Hi,h2,G) in 



12) will be 

maximized by choosing G to be a diagonal matrix, i.e. G = A G ; 
that is equivalent to choosing U G = I . Considering that the spectral 
decomposition of G and G in 



13 1 are G = 



UgA-gUg and 
J in {dI 



U G = 



G = U qKqU g , respectively, then choosing 
dictates that U g = Us must hold. 

So far, it is proved that A G can achieve the same capacity as 
G. However, for a complete proof, it must be made sure that by 



choosing G = A G the power constraint of the relay defined in 1 14 1 
is not violated. It is clear from l |14| > that the power constraint of the 
relay depends on G only via two terms: T^A^G) and Tr(A E 2 G). 
In what follows, we will prove that by choosing G = A G , both terms 
in the relay power constraint will be minimized, hence, fulfilling the 
constraint. 

Lemma: Let G be an arbitrary positive symmetric matrix with 
spectral decomposition G = U G A G U^,. Let As be a diagonal 
matrix with the elements on the diagonal in descending order and 
k > 0. Then 



Tr(A 



£ A G. 



< Tr(A E fc G) 



(17) 



Proof: See (17| Theorems 1-3] for a detailed proof for k — 1. 
The proof can easily be extended to arbitrary k > 0. ■ 
By applying \17\ to ( [14) , the power constraint will be fulfilled; that 
is, by choosing Ug = Us, or equivalently G = A G — AqA|, it 
is assured that 



PsTr(A^A 6 ) + iV Tr(A E 2 A 6 ) < 



P s Tr(A E 1 G) + iVoTi(Aj G) < P R . 



(18) 



By completing the proof for the power constraint of the relay, the 
answer to Problem 1 in Section |n] is accomplished. Therefore, the 
transmission from the relay should be conducted in the direction of 
the eigenvectors of the correlation matrix E. However, the power per 
eigenvector (Xf) must be determined using numerical methods (see 
e.g. 1 1 8| for a single hop system model). 

So far, we assumed equal correlation matrix E in Hi and hi 
which implies that the correlation is a consequence of space limit 
in the relay node, however, when the correlation occurs due to 
unobstructed relay node, the correlation matrices corresponding to 
Hi and h,2 can be different. Nevertheless, derivation of F will not be 
much different. Let assume that Hi = Ej Hi w and h,2 = /i2™S|. 
Then, applying a similar method as explained in this section, one can 
prove that the optimal relaying matrix is 

tH 



F = Us 2 AU, 



(19) 



where Us i , i £ {1, 2}, is a unitary matrix corresponding to the 
eigenvectors of Ej and A is a diagonal matrix with its elements in 
decreasing order which is to be calculated numerically. 

IV. Optimality of Maximum Eigenmode Relaying 



In this Section, Problem 2 from Section |II-C| is addressed. The 
focus of this section is to derive necessary and sufficient conditions 
under which MER is the optimal transmission method from the relay. 

By the results obtained in Section [ill] and considering that A Q = 



AgA e , the ergodic capacity in {7) can be rewritten as 
1 



E{C(Hi w ,h 2w ,A G )} (20) 



C — - max 

A G :P s Tr(A s A G )+iV Tr(A G )<P R 

where the expectation-operation is carried out over Hi w and h 2n 
with 

C(Hi w , h 2w ,A G ) = log (1 + 



^h 2w A G As H i w Hf w A a Ash2 W , 
n s (l + h 2w A G A'ih^ m ) 
Then C(Hi w , h,2w, Ag) can be simplified according to 



(21) 



7 I h,2 W A G A s hi v 



C(H lw ,h 2w , A G ) = log (1 + 



ns(l + h 2w A G A s h2 W ) 



where hi W: i represents the ith column of Hi w . One can write the 
numerator 



1 " s . 1 nR 

- ]T I h 2w AlAshi w , l | 2 = ]T Xf\f 2 XjY 



ns . 

1—1 

and the denominator 



(23) 



3=1 



h 2w A G Alh? w = Y, ^f x T X J (24) 

3=1 

where Xj =\ h2w,j j 2 is an exponential random variable with unit 
mean, i.e. fxj{t) = e~ t , and Y — X^i=i ! \ 2 is tne sum 

of ns i.i.d exponential random variables with parameter ns. Indeed, 
Y has an Erlang-distribution with rate and shape equal to ns, i.e. 
f Y (t) = {t ns - l)e~7(n s - 1)!. The proof for the identity |23}is 
omitted due to lack of space but a similar proof can be found in 1 12 
Eq. (71)-(72)]. Plugging in |23j and <[24]> into {22), it will simplify 
to 



C(Hi w ,h 2w , A G ) = log (l 



1 + 



(25) 



"R "R 

= log + Y Af Af (1 + iY)Xj) - (1 + E X ? X T X i) ■ 

3=1 3=1 



A. Necessary Condition 



Now we turn the attention to the question under which condition 
MER is optimal; i.e. the condition under which = for j > 2 is 
the capacity achieving transmission method with 

Pr 



X? 



(i + a?p s ; 



(26) 



which is calculated by substituting Xj — for j > 2 in 

In order to evaluate the necessary and sufficient condition on the 
optimality of MER, we assume P = 2~2j=i ■ Assume that the 
power P—p is allocated to the dominant eigenvector of G, i.e. = 
P — p, and the power p is allocated to the remaining eigenvectors 
of G, i.e. p = J2712 A? • H is clear that if MER is optimal, then 
dC'(p)/dp |p=o< which provides the necessary condition for the 
optimality of MER and d 2 C(p)/dp 2 |p=o< which confirms the 
sufficiency. It can be proved that dC'(p)/dp \ p= o will be maximized 
if X 2 — p and Xf = for j > 2; similar discussions are provided 
in 1 1 2| . |17| , |19| and, hence, we do not repeat it here. Therefore, 
given the assumption Xf — P — p and A^ = p, \25\ will be further 
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simplified to 

C(H lw ,h 2w ,A G )= (27) 

log (1 + (P - p)aP(i + 7 y)Xi + paP(i + 7 y)x 2 



- log I 1 + (P - p)Af Xx + pAf 2 X 2 



and so, dC(p)/dp \ p= o from 1 20 1 and {27} equals 

AP(i + 7 y)x 2 



9p 



E{C(iTi M) fo 2K „A G )} | p= o = E{ 



i + PAp(i + 7 y)X! 
AP(i + 7 y)^! 



} 



1 1 + PAP(1 + 7 F)X/ 



The first expectation in {28} can be written as 



E{ 



AP(i + 7 y)x 2 



} 



AfE{ : 



(i + 7 y) 



} (29) 



i + PAp(i + 7 y)xP 2 L i + PAp(i + 7 y)Xr 

which follows from the fact that X 2 is independent of Xi and Y and 
that E{X 2 } = 1. By simple manipulations, the second expectation 
in \2S\ equals 



E{ 



AP(i + 7 y)Xx 
i + PAP(i + 7 y)Xi 
i 

p 



} = 



P i l + PAP(l + 7 y)X 1 



and the last expectation in {28) can be written as 



E{- 



APXi - 



1 + PAPXi 



(30) 



} 



(3D 



(i + PAP)r(o,- 



+ 



P 2 AP 

where r(0, z) is the incomplete Gamma function |2()| Sec. 6.5]. 
By combining {20}, ( |28) , \29\ , l |30) and < | 3 X | > , in order to compute 

dC(p)/dp | p =o< 0, we have 



dp 



\p=o 



E{ 



AP(i + 7 y) + i/P 
i + PAp(i + 7 y)Xi 



} 



(32) 



(i + PAP)r(o,^ 



< 



(PAP 2 

from which by elementary operations the following constraint will 
be obtained for the MER to be optimal: 



aP< 



! r(o, 



i_ 



-paPe{^} 



(33) 



s2Tcr(i+7y) 



P 2 APE{ 



} -Pe J 



T(0, 



where Z=l + P\f 2 (l + jY)Xl Note that E{l/Z} and E{(1 + 
7 y ) /Z} in \33\ are non-trivial and do not seem to have closed-form 
solutions. However, they can be calculated either using monte carlo 
simulations or by numerical integration of 



with 



E 4> = 



A(t) 



A(t)e A(t) T(0,A(t))dt 



(34) 



PAf 



T(0,A(t))dt (35) 



(36) 



PAP(l + 7 i) ' 

The proof of |34} and \35\ follows from the definition of the 



expectation operation and is omitted due to lack of space. 

B. Sufficient Condition 

By deriving the necessary condition for MER to be optimal, it can 
be proved that the necessary condition is also a sufficient condition for 
MER to be optimal. It will be proved by showing that d 2 C(p) / dp 2 < 
for arbitrary p £ [0, P] 



For deriving d 2 C(p)/dp 2 , from |20| and 1 27 we have 

dp 2 
Ps((A 



K{C(H lw ,h 2w ,A a )} \ P=0 = (37) 

x 1 - (Aj) 2 x 2 ) 2 (2 + 7 y + 2P(AP 2 Xi(i + 7 y))y 



vP 2 



nsiV (l + P(Ap2Xx) 2 (l + P(Ap 2 Xx(l +jY)) z 

It is clear from (37} that d 2 C{p)/dp 2 < for every p, and so, the 
necessary condition for the optimality of MER derived in {33} is 
sufficient as well. 



V. Discussion 

In the previous section, a necessary and sufficient condition were 
derived which show the optimality region of MER. As it is clear 
from |32j or (33) , there is need to perform Monte Carlo simulations 
or numerical integrations to validate the optimality of MER. However, 
computationally expensive Monte Carlo simulations and also numer- 
ical integration methods are not a good choice in practical real-time 
communication systems; therefore, in the sequel, we will investigate 
two simplified approaches which lead to closed form solutions. 

We first derive a lower bound for the necessary and the sufficient 
conditions for the optimality of MER which are a direct result 
of lensen's inequality. Consequently, if the lower bound condition 
confirms the optimality of MER, there is no need for the Monte Carlo 
simulations and, hence, MER can be used as a capacity-achieving 
method. 

Then we investigate the effect of the source antenna array on the 
optimality of MER. That is done by considering large antenna arrays 
in the source using the central limit theorem. Novel, closed form 
necessary and sufficient conditions are derived for the optimality of 
MER which, according to the simulations, turn out to very tightly 
approximate the results for arbitrary numbers of antennas in the 
source. 



A. A Lower Bound on the Optimality of Maximum Eigenmode 
Relaying 

From previous discussions, it is clear that dC(p)/dp < | p= o 
defines the condition under which MER is optimal. We will exploit 
it to derive necessary and sufficient conditions which specify a lower 
bound on the optimality of MER. From careful inspection of {32} it 
is clear that 

AP(i + 7 y) + i/p 



> 



(38) 



i + PAp(i + 7 y)Xi 

which equals the expression inside the expectation operation in 
(32) . On the other hand, according to lensen's inequality whe have 
E{ij(X)} > gr(E{X}) for a convex function g(x). Let us assume 
g(X) — 1/X, which is a convex function for X > 0. Therefore, by 
applying Jensen's inequality to the expectation operation in {32}, we 
have 



E{ 



AP(i + 7 y) + i/p 



1 + PA? 2 (1+7*0*1 



}> 



(39) 



l + PAPq+^Xg 

1 1 Ap(i + 7 y) + i/p ] 
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Fig. 1. MER optimality region for (Ps.Pr) pairs for various correlations S, 



with 

i + PAP(i + 7 m 1 

1 AP(l + 7V) + l/^ 



(40) 



P(AP(1 + PAP) + n s (AP - AP)e p E„ s+1 (P)) 
AP(1 + PAP) 

where 2? = n s (l + PAP)/7PAP and E m (z) the Exponential 
Integral function [20 Sec. 5.1]. Therefore^] the following expression 
is defined as a lower bound which specifies that MER is optimal: 



paPaP(i + pa| 



AP(1 + PAP) + n s (AP - AP)e B E rls+1 (D) 



< 



(41) 



e PA r (i + pAP)r(o, 



PAf 



B. Large Antenna Array in the Source 

It is interesting to investigate the effect of the dimension of the 
antenna array on the optimality of MER. Note that as we focus on 
the optimality of maximum eigenmode transmission from the relay; 
therefore we only consider the effect of large antenna arrays on 
the optimality of MER. Assuming ns 3> 1, the random variable 
Y ~ E(ns, ns) can be approximated by a Gaussian distribution with 
mean 1 and variance 1/ns, i.e., approximately, Y ~ 7V(l,l/ns). 
Therefore, as ns — > oo, we obtain Y — ¥ 1. By substituting Y = 1 
in {33}, we have 1 + 7Y = 1 + 7. Hence, E{l/Z} in {33) can be 
written in closed form as 



E{-}=Aie Al F(0,Ai) 



(42) 



where Ai = l/PAp(l + 7) is obtained by substituting t = 1 in 
l |36[ l. Then, the necessary and sufficient conditions for the optimality 
of MER, i.e. dC(p)/dp | p=0 < in {32) or {33) for large n s , after 



some manipulation lead to the following closed-form constraint: 



aP< 



T(0, 



P\f 2 Aie Al T(0,Ai)} 



(43) 



Pe A ^T(0,Ai) -Pe' 



! T(0, 



'The proof of the expectation operation in |40J is omitted due to space 
limit, however, it can be validated using symbolic tools, e.g. Mathematica. 



Fig. 2. Mutual information between input and output of the relay system 
with rig = 2 = ?ir and MER. Given Pg = lOdB, the Pr corresponding to 
solid lines illustrates the region where MER achieves capacity, while dashed 
lines show the region where MER can not achieve capacity. 



It is clear from §43\ that with large antenna array in the source, the 
optimality of MER depends on the source transmission only via Ps 
and it is independent of H\, 

It will be shown, by numerical simulations in the next Section, 
that the closed-form constraint in {43) approximates {33) with high 
accuracy. 

VI. Numerical Results 

Computer simulations are provided to verify the analytical results 
derived in the previous sections. We assume that the relay is equipped 
with two antennas, but various numbers of the antennas in the source 
are evaluated. As explained in Section [TTJ we assume i.i.d Hi ]W 
and h,2. w , where the elements are circularly symmetric complex 
Gaussian random variables with zero mean and unit variance (block 
Rayleigh fading assumption). The noise power in the relay and the 
destination is assumed to be A^o = 1. Fig. [T| illustrates the (Ps, Pr)- 
pairs for which MER is optimal. The optimality region of MER is 
calculated for p — 0.3 and 0.5 where p is defined as the inter- 
antenna correlation or the off-diagonal elements of the correlation 
matrix S. It is clear from Fig. [T] that by increasing p, the optimality 
region of MER increases, too. Another interesting observation from 
Fig. [T] is that the optimality of MER is almost independent of the 
number of antennas in the source node. The figure illustrates the 
optimality region using the expression derived in {33) for ns = 2 
and 4 which require numerical integrations and also using the closed 
form expression derived in {43| when ns — ► 00; it explicitly shows 
that regardless of ns, the optimality regions of MER coincide with 
very little difference. Fig. [2] shows the mutual information between 
the input and output of the specified MIMO relay channel. Note that 
the solid lines achieve capacity using MER but the dashed lines do 
not achieve capacity and, hence, MER is not optimal in this part of 
the curve. 

VII. Conclusion 

A dual hop cooperative system was investigated in this paper. The 
source and the relay nodes are equipped with multiple antennas and 
the destination with single antenna. Optimal precoding matrix in the 
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relay was derived. It was shown that the optimal transmission from 
the relay should be conducted in the direction of the eigenvectors 
of the transmit-channel covariance matrix. Then, a necessary and 
sufficient condition was derived, under which, the relay transmission 
only from the strongest eigenvector achieves capacity; this method of 
transmission was called Maximum Eigenmode Relaying. The exact 
result contains two integrals which need to be solved numerically. 
Moreover, a closed form lower bound was derived for the optimality 
region of MER. We further investigate to evaluate the effect of the 
source antenna array on the optimality of MER and derived closed 
form expression when the source is equipped with infinite antennas. 
The simulation results show that MER optimal region with infinite 
antennas in the source coincides with MER optimal region when 
the source is equipped with much lower number of antennas with 
inappreciable difference. 
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